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1. INTRODUCTION 

At present several actual problems of the solid-state electronics (such as increasing performance, 
reliability, and density of elements of integrated circuits: diodes, field-effect, and bipolar transistors) are 
intensively solving [1]-[6]. To increase the performance of these devices, it is attracted an interest 
determination of materials with higher values of charge carrier's mobility [7]-[10]. One way to decrease the 
dimensions of elements of integrated circuits is by manufacturing them in thin-film heterostructures [3]-[6]. 
In this case, it is possible to use inhomogeneity of heterostructure and necessary optimization of doping of 
electronic materials [11] and development of epitaxial technology to improve these materials (including 
analysis of mismatch induced stress) [12]-[14]. An alternative approach to increasing the dimensions of 
integrated circuits is the use of laser or microwave types of annealing [15]-[18]. 

The main objective of the present paper is to formulate a condition to decrease the dimensions of 
elements (first of all field-effect heterotransistors) of a nulling-resistor output amplifier as shown in Figure 1, 
Figure 1(a) shows the structure of the considered nulling-resistor output amplifier and Figure 1(b) shows 
heterostructure with a substrate, epitaxial layers and buffer layer (view from the side), manufactured in the 
framework of a heterostructure with a spatial configuration: we consider a heterostructure, which consists of 
a substrate and an epitaxial layer; the epitaxial layer includes into itself several sections, which were 
manufactured by using other materials; the obtained sections were doped by diffusion or by ion implantation 
to obtain the required type of conductivity (p or n) and used as sources, drains and gates of field-effect 
transistors. The accompanying aim of the present paper is an analysis of the possibility to increase the density 
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of these transistors. Heterostructures have mismatch-induced stress due to the difference between values of 
lattice distances of the substrate and the epitaxial layer. The stress could lead to dislocation inconsistencies in 
the substrate and epitaxial layer. In the framework of the paper, we formulate conditions to decrease 
mismatch-induced stress. The decrease is another accompanying aim of the present paper. To make a 
prognosis of manufacturing of the considered amplifier we introduce an analytical approach for the prognosis 
of mass and heat transport during technological processes without crosslinking of solutions on interfaces 
between layers in the considered heterostructure. The introduced analytical approach gives a possibility to 
take into account the dependence of parameters of technological parameters on spatial coordinates as well as 
the dependence of these parameters on time and at the same time nonlinearity of considered processes. 
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Figure 1. Formulation of condition to decrease dimensions of the elements of (a) Structure of the considered 
nulling-resistor output amplifier and (b) Heterostructure with a substrate, epitaxial layers, and buffer layer 
(view from the side) 


2. METHOD OF SOLUTION 

To solve the main aim of the present paper we consider a heterostructure, which consists of a 
substrate and an epitaxial layer as shown in Figure 1. We also consider a buffer layer between the substrate 
and the epitaxial layer. The epitaxial layer includes into itself several sections, which were manufactured 
using another material. These sections have been doped by diffusion or ion implantation to manufacture the 
required types of conductivity (p or n). These areas became sources, drains, and gates as shown in Figure 1. 
After the doping, is required annealing of dopant and/or radiation defects. We analyzed the redistribution of 
dopant and radiation defects to determine conditions, which correspond to decreasing dimensions the 
considered of the field-effect transistors in the framework of the considered amplifier and at the same time 
increasing their density. At the same time, we consider the possibility to decrease mismatch-induced stress. 
To make the analysis we determine and analyze the distribution of concentration of dopant in space and time 
in the considered heterostructure. We determine the distribution by solving the second Fick's law in the 
following form [1], [19]-[23] with boundary and initial conditions. 
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Where C (x,y,z,f) is the Spatio-temporal distribution of concentration of dopant; © is the atomic 
volume of dopant; V; is the symbol of surficial gradient; I C(x, y, z, t)dz is the surficial concentration of 
dopant on the interface between layers of heterostructure (in this situation we assume, that Z-axis is 
perpendicular to the interface between layers of heterostructure); zu(x,y,z,t) is the chemical potential due to 
the presence of mismatch-induced stress; D and Ds are the coefficients of volumetric and surficial diffusions. 
Values of dopant diffusions coefficients depend on properties of materials of the heterostructure, speed of 
heating and cooling of materials during annealing, and Spatio-temporal distribution of concentration of 


dopant. Dependences of dopant diffusions coefficients on parameters could be approximated by the following 
relations [21]-[23]. 
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Where Dz (x,y,z,T) and Dzs (x,y,z,T) are the spatial (due to accounting for all layers of heterostructures) and 
temperature (due to Arrhenius law) dependences of dopant diffusion coefficients; T is the temperature of 
annealing; P (x,y,z,7) is the limit of solubility of dopant; parameter y depends on properties of materials and 
could be an integer in the following interval » e[1], [3], [21]; V (x.y,z.t) is the Spatio-temporal distribution of 
concentration of radiation vacancies; V' is the equilibrium distribution of vacancies. The concentration 
dependence of the dopant diffusion coefficient has been described in detail in [21]. Spatio-temporal 
distributions of concentration of point radiation defects have been determined by solving the following 
system of equations [19], [22], [23] with boundary and initial conditions. 
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Where J (x,y,z,f) is the spatio-temporal distribution of concentration of radiation interstitials; J" is the 
equilibrium distribution of interstitials; Di(x.y,z, T), Dv(x,y,z,T), Dis(x.y,z,T), Dvs(x,y,z,T) are the coefficients 
of volumetric and surficial diffusions of interstitials and vacancies, respectively; terms V7(x,y,z,f) and 
P(x,y,z,t) correspond to generation of divacancies and diinterstitials, respectively (see, for example, [23] and 
appropriate references in this book); krv(x,y,z,T), kiAx,y,z,T) and kv,v(x,y,z,T) are the parameters of 
recombination of point radiation defects and generation of their complexes. Spatio-temporal distributions of 
divacancies ®,(x,y,z,t) and diinterstitials ®(x,y,z,t) could be determined by solving the following system of 
equations [19], [22]—-[24]. 
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With boundary and initial conditions, 
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Where Da(x,y,z,T), Dov(x,y,z,T), Dars (x,y,z, T) and Doys(x,y,z,T) are the coefficients of volumetric 
and surficial diffusions of complexes of radiation defects; k)(x,y,z,T) and kv(x,y,z,T) are the parameters of 
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decay of complexes of radiation defects. Chemical potential 4; in (1) could be determine by the following 
relation [19]. 
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Where E(z) is the Young modulus, oj is the stress tensor; uj; = (2 + u) is the deformation tensor; ui, uj 
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are the components ux(x,y,z,f), Uy(x,y,z,t) and u-(x,y,z,f) of the displacement vector U(x, y, zZ, t); xi, x; are the 
coordinate x, y, z. The (3) could be transformed to the following form, 
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Where o is the Poisson coefficient; & = (a;-acr)/agr is the mismatch parameter; ds, az; are lattice distances of 
the substrate and the epitaxial layer; K is the modulus of uniform compression; £ is the coefficient of thermal 
expansion; 7. is the equilibrium temperature, which coincides (for our case) with room temperature. 
Components of displacement vector could be obtained by solution of the following equations [25]. 
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x B(Z)[T(~ y,z, t) — T,], p (z) is the density of materials of the heterostructure, A: Is the Kronecker symbol. 
With an account, the relation for oj; last system of the equation could be written as, 
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Conditions for the system of (8) could be written in the form. 
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We determine Spatio-temporal distributions of concentrations of dopant and radiation defects by 
solving the (1), (3), and (5) framework standard method of averaging function corrections [26]. Previously 
we transform the (1), (3), and (5) to the following form with an account of initial distributions of the 
considered concentrations, 
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Farther we replace concentrations of dopant and radiation defects on the right sides of (1a), (3a), and 
(5a) on their not yet known average values ap. In this situation, we obtain equations for the first-order 
approximations of the required concentrations in the following form, 
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Integration of the left and right sides of the (1b), (3b), and (5b) on time gives us the possibility to 
obtain relations for approximation in the final form 
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We determine average values of the first-order approximations of concentrations of dopant and 
radiation defects by the following standard relation [26]. 
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Substitution of the relations (1c), (3c), and (5c) into relation (9) gives us the possibility to obtain the 
required average values in the following form 
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Where we determine approximations of the second and higher orders of concentrations of dopant and 
radiation defects framework standard iterative procedure of a method of averaging of function corrections 
[26]. Framework this procedure to determine approximations of the n-th order of concentrations of dopant 
and radiation defects we replace the required concentrations in the (1c), (3c), (5c) on the following sum 
Onp*p n-1(xX,y,2,t). The replacement leads to the following transformation of the appropriate equations 
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Integration of the left and the right sides of (1d), (3d), and (5d) gives us the possibility to obtain 
relations for the required concentrations in the final form 
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Average values of the second-order approximations of required approximations by using the 
following standard relation [26]. 
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Substitution of the relations (1e), (3e), and (5e) into relation (10) gives us the possibility to obtain 
relations for required average values a 2p. 
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Farther we determine solutions of (8), i.e. components of the displacement vector. To determine the 
first-order approximations of the considered components framework method of averaging function 
corrections we replace the required functions on the right sides of the equations with their not yet known 
average values ai. The substitution leads to the following result. 
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Integration of the left and the right sides of the relations on-time t leads to the following result, 
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Approximations of the second and higher orders of components of the displacement vector could be 
determined by using standard replacement of the required components on the following sums ai+ui(x,y,z,t) 
[26]. The replacement leads to the following result. 
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Integration of the left and right sides of the relations on-time t leads to the following result, 
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In the framework of the paper, we determine the concentration of dopant, concentrations of radiation 
defects, and components of displacement vector by using the second-order approximation framework method 
of averaging function corrections. This approximation is usually enough good approximation to make a 
qualitative analysis and to obtain some quantitative results. All obtained results have been checked by 
comparison with the results of numerical simulations. 


3. RESULTS AND DISCUSSION 

In this section, we analyzed the dynamics of redistributions of dopant and radiation defects during 
annealing and under influence of mismatch-induced stress. Typical distributions of concentrations of dopant 
in heterostructures are presented in Figures 2 and 3 for diffusion and ion types of doping, respectively. These 
distributions have been calculated for the case when the value of the dopant diffusion coefficient in the 
epitaxial layer is larger than in the substrate. The figures show, that inhomogeneity of heterostructure gives 
us the possibility to increase the sharpness of p-n- junctions. At the same time, one can find increasing 
homogeneity of dopant distribution in the doped part of the epitaxial layer. Increasing of sharpness of the p- 
n-junction gives us the possibility to decrease switching time. 

The second effect leads to decreasing local heating of materials during the functioning of the p-n- 
junction or decreasing of dimensions of the p-n-junction for a fixed maximal value of local overheat. 
However, in the framework of this approach to the manufacturing of bipolar transistors, it is necessary to 
optimize the annealing of dopant and/or radiation defects. The reason for this optimization is the following. If 
the annealing time is small, the dopant did not achieve any interfaces between materials of the 
heterostructure. In this situation, one cannot find any modifications in the distribution of concentration of 
dopant. If the annealing time is large, the distribution of concentration of dopant is too homogenous. We 
optimize the annealing time framework recently introduces approach [15], [25]-[32]. To framework this 
criterion we approximate the real distribution of concentration of dopant by step-wise function as shown in 
Figure 4 and Figure 5. Farther we determine optimal values of annealing time by minimization of the 
following mean-squared error. 

Where y (x, y, z) is the approximation function. Dependences of optimal values of annealing time 
on parameters are presented in Figures 6 and 7 for diffusion and ion types of doping, respectively. It should 
be noted, that it is necessary to anneal radiation defects after ion implantation. One could find the spreading 
of the concentration of distribution of dopant during this annealing. In the ideal case distribution of dopant 
achieves appropriate interfaces between materials of heterostructure during the annealing of radiation defects. 
If the dopant did not achieve any interfaces during annealing of radiation defects, it is practicably to 
additionally anneal the dopant. In this situation, the optimal value of the additional annealing time of 
implanted dopant is smaller, than the annealing time of infused dopant. 

Farther we analyzed the influence of relaxation of mechanical stress on the distribution of dopants in 
doped areas of the heterostructure. Under the following condition &< 0 one can find compression of 
distribution of concentration of dopant near the interface between materials of the heterostructure. Contrary 
(at &>0) one can find a spreading of distribution of concentration of dopant in this area. This change of 
distribution of concentration of dopant could be at least partially compensated by using laser annealing [29]. 
This type of annealing gives us the possibility to accelerate the diffusion of dopants and other processes in 
the annealed area due to the inhomogenous distribution of temperature and Arrhenius law. Accounting 
relaxation of mismatch-induced stress in heterostructure could lead to changing of optimal values of 
annealing time. Mismatch-induced stress could be used to increase the density of elements of integrated 
circuits. On the other hand, could lead to generation dislocations of the discrepancy. Figure 8 shows 
distributions of components of the displacement vector, which is perpendicular to the interface between 
layers of the heterostructure. 
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Which is perpendicular to the interface between the epitaxial layer substrate. Increasing the number 
of curves corresponds to increasing of difference between values of the dopant diffusion coefficient in layers 
of heterostructure under conditions when the value of the dopant diffusion coefficient in the epitaxial layer is 
larger than the value of the dopant diffusion coefficient in the substrate. 
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Figure 2. Distributions of concentration of infused dopant in heterostructure from Figure 1 in direction 
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Figure 3. Distributions of concentration of implanted dopant in heterostructure from Figure 1 in direction 


Which is perpendicular to the interface between the epitaxial layer substrate. Curves 1 and 3 
corresponds to annealing time © = 0.0048(Lx2+Ly2+Lz2)/DO0. Curves 2 and 4 corresponds to annealing time 


© = 0.0057(Lx2+Ly2+Lz2)/D0. Curves 1 and 2 correspond to a homogenous sample. Curves 3 and 4 
correspond to heterostructure under the condition when the value of dopant diffusion coefficient in the 


epitaxial layer is larger than the value of dopant diffusion coefficient in the substrate. 
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Figure 4. Spatial distributions of dopant in heterostructure after dopant infusion 
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Curve 1 is the idealized distribution of dopant. Curves 2-4 are real distributions of dopants for 


different values of annealing time. Increasing the number of curves corresponds to increasing of annealing 
time. 
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x 
Figure 5. Spatial distributions of dopant in heterostructure after ion implantation 


Curve 1 is the idealized distribution of dopant. Curves 2-4 are real distributions of dopants for 


different values of annealing time. Increasing the number of curves corresponds to an increase in annealing 
time. 
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Figure 6. Dependences of dimensionless optimal annealing time for doping by diffusion 


Which have been obtained by minimization of mean-squared error, on several parameters. Curve 1 
is the dependence of dimensionless optimal annealing time on the relation a/L and € = y = 0 for equal to each 
other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of 
dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and € = y = 0. Curve 3 is the 
dependence of dimensionless optimal annealing time on value of parameter € for a/L=1/2 and e = y = 0. 


Curve 4 is the dependence of dimensionless optimal annealing time on value of parameter y for a/L=1/2 and € 
= & = 0. 
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Figure 7. Dependences of dimensionless optimal annealing time for doping by ion implantation 


Which have been obtained by minimization of mean-squared error, on several parameters. Curve 1 
is the dependence of dimensionless optimal annealing time on the relation a/L and € = y = 0 for equal to each 
other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of 
dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and € = y = 0. Curve 3 is the 
dependence of dimensionless optimal annealing time on value of parameter € for a/L=1/2 and e = y = 0. 
Curve 4 is the dependence of dimensionless optimal annealing time on value of parameter y for a/L=1/2 and € 
-6z0. 


Figure 8. Normalized dependences of component uz of displacement vector on coordinate z for nonporous 
(curve 1) and porous (curve 2) epitaxial layers 


4. CONCLUSION 

In this paper, we model redistribution of infused and implanted dopants with account relaxation 
mismatch-induced stress during manufacturing field-effect heterotransistors in the framework of a nulling- 
resistor output amplifier. We formulate recommendations for optimization of annealing to decrease 
dimensions of the considered transistors and to increasing of their density. In this paper, we formulate a 
recommendation to decrease mismatch-induced stress by using radiation processing of stressed areas of the 
heterostructure. This effect was recently found experimentally, but in this paper, we introduce a model of this 
effect and an analytical approach to the prognosis of this effect. We also introduce an analytical approach for 
prognosis of diffusion and ion types of doping with account concurrent changing of parameters in space and 
time has. At the same time, the approach gives us the possibility to take into account the nonlinearity of 
considered processes. The considered approach is the significantly more common approach in comparison 
with other ones for related applications. At the same time, analytical approaches to the solution of boundary 
problems give explicit functional dependences of obtained solutions on all parameters in comparison with 
numerical ones. 
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